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Classically, the Prigogine-Defay (PD) ratio involves differences in isobaric volumic 
specific heat, isothermal compressibility and isobaric thermal expansion coefficient 
between a super-cooled liquid and the corresponding glass at the glass transition. 
However, determining such differences by extrapolation of coefficients that have been 
measured for super-cooled liquid and glassy state, respectively, poses the problem that 
it does not take into account the non-equilibrium character of the glass transition. 
In this paper, we asses this old question by taking into account the gradual change of 
configurational contributions to the three thermodynamic coefficients upon varying 
temperature and pressure. Macroscopic non-equilibrium thermodynamics is applied 
to obtain a generalized form of the PD ratio. The classical PD ratio can then be taken 
as a particular case of this generalization. Under some assumptions, a configurational 
PD ratio (CPD ratio) can be expressed in terms of fictive temperature and fictive 
pressure what, hence, provides the possibility to experimentally verify this formalism. 
Noteworthy and differing from previous approaches towards the PD ratio, here, the 
glass transition is considered as non-isoaffine. 
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I. INTRODUCTION 

Upon vitrification, some of tlie tliermodynamic properties of a glass former undergo a jump 
which can be observed experimentally. Classically, the Prigogine-Defay (PD) ratio was defined as 
the ratio of the jumps of the following thermodynamic quantities: 

where ACp, An^ and Aa^ are the differences taken at the glass temperature, Tg , of the extrapol- 
ated measured isobaric heat capacity, isothermal compressibility and isobaric thermal expansion 
coefficient as observed in liquid and glassy states, respectively. V is the volume of the system. 
The PD ratio has initially been derived by Prigogine and Defay, using classical thermodynamics of 
irreversible processes and De Bonder's order parameter framework l|. In their original definition, 
the sign "A" represents a difference between a state at constant affinity (isoaffine state) and a 
state at constant order parameter (that is, the glassy state). Indeed, since the affinity and the 
order parameter are thermodynamically conjugated, in analogy to the well-known expression of the 
difference between isobaric and isochoric heat capacity {Cp — Cy), an expression of the difference 
between isoaffine and iso-order-parameter heat capacity {Ca — C^) can be derived. So, in allowing 
the affinity to tend towards zero, the sign "A" represents the super-cooled liquid-to-glassy differ- 
ence and, hence, the common definition of the PD ratio is obtained. If applied to second order 
phase transitions, simple equilibrium thermodynamic arguments show that this quantity is equal 
to unity. This can be demonstrated by expressing the continuity of volume and entropy (first 
order derivatives of the Gibbs free energy) following Ehrenfest's classification, when the system 
undergoes a phase change 

For the case of glass transitions, however, numerous experiments have revealed that this ratio 
may differ from unity. That is, in some cases, PD < 1 has been observed, but more often, 
experimental value of PD > 1 are obtained (as an extreme example, for the case of vitreous silica, 
classical experimental approaches yield a value of PD > 10000). In a pioneering study, Davies and 
Jones were transferring De Donder's nomenclature to the classical field of glass sciencejsl. They 
interpreted the experimental observation of PDt^I as an evidence that multiple (more than one) 
order parameters must be introduced to completely describe the glass transition in a real system. 
Over the following decades, the so-called "order parameter approach" has then been applied and 
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refined by Moynihan, DiMarzio, Goldstein, Gupta and others 4t8| and became widely accepted. 
A more general approach towards the PD ratio was undertaken by Lesikar and Moynihan. In their 
consideration, the value of the PD ratio results from a thermodynamic stability condition, and can 
be generalized to unconventional thermodynamic parameters such as electric polarisation ^\. It 
was assumed that each thermodynamic variable is governed by its corresponding order parameter, 
what generally implies PD > 1. They also have derived a particular condition for which the PD 
ratio is equal to unity even if several order parameters are taking part in the glass transition 9|. 



The topic has been reviewed in a very complete form by Nemilov |10|. In further studies, also, 
Gutzow and Schmelzer have discussed the PD ratio mostly from viewpoint of its genuine definition 

HI- 

During the last decade, however, with the publication of several new papers, the debate has 
been reopened. Similar to Ritland's original concept of the fictive temperature [l^, a phe- 



nomenological "two temperature thermodynamic approach" was developed by Nieuwenhuizen in 
which he reconsiders various classical problems of the glass transition 1^, [l5|. Among them is 



counted the non-unity of the PD ratio. It was shown that by taking into account the config- 
urational entropy of a supercooled system, one of the two Ehrenfest relations is modified what, 
consequently yields PD ^ 1. The problem was approach from another angle by Ellegaard and 
al. [l6|. In their study, they are introducing a dynamic PD ratio where only the imaginary parts 
of the involved complex thermoviscoelastic coefficients are taken into account. They demonstrate 
rigorously that under such assumptions PD may become unity at all frequencies even when only a 
single order parameter is taken into account. Experimental and theoretical works have confirmed 



this for volume-enthalpy correlated glass-formers 
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18|. Eventually, in 2006, Schmelzer and 



Gutzow, in re-formulating a thermodynamic description by means of classical thermodynamics of 



irreversible processes, derived a new expression of the PD ratio jl9|. A new expression of the P.D. 
ratio was derived. Their expression highlights the role of thermodynamic affinity, explicitly stating 
that the considered system is out of thermodynamic equilibrium, and that the glass transition 
does not follow Ehrenfest's original definition of a phase transition (where a discontinuity in some 
derivative of Gibbs free energy occurs at a precise temperature or pressure). More precisely, the 
transition process takes place over a definite temperature or pressure range. Within this range, 
the system starts to deviate from thermodynamic equilibrium. Thus, in their formulation, the 
difference "A" attains its genuine formulation: a difference between a state at constant affinity 
and a state at constant order parameter. In the present study, based on these considerations, a 
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simple expression of the PD ratio is derived, based solely on enthalpy and affinity. Subsequently, 
it is demonstrated within this formulation, that PD can differ from unity, independently on the 
number of involved order parameters. 

Under the "single order parameter assumption", we provide a generalization of Schmelzer and 
Gutzow's expression of the PD ratio, and we extend the approach to other than isoaffine transform- 
ations. The configurational parts of the three thermodynamic coefficients, Cp, kt and ap are used 
directly in order to define a configurational Prigogine-Defay (CPD) ratio. The classical expression 
of the PD ratio as well as the expression of Schmelzer and Gutzow become particular cases of 
this definition, where the creation of affinity occurs as a direct result of changes in temperature 
or pressure. This demonstrates the complexity of the phenomenological behaviour of the glass 
transition as a thermodynamically irreversible process. 

The paper is written as follows: In section 2, the non-equilibrium thermodynamic approach of 
the glass transition is briefly exposed. In section 3, the CPD ratio is proposed. A tri-dimensional 
diagram is then depicted and the glass transition is seen as a path in this diagram. A simple 
expression of the CPD ratio is then provided, revealing only the isobaric heat capacity and the 
isothermal compressibility coefficient. Then, this last expression is simply transformed to make 
apparent the fictive temperature and the fictive pressure. In section 4, the unity of the CPD ratio 
is discussed. 



II. NON-EQUILIBRIUM THERMODYNAMIC APPROACH OF THE GLASS 

TRANSITION 



A. The order parameter model 

A physical system at thermodynamic equilibrium is characterized by variables or parameters 
defining the thermodynamic state of the system. These variables {xj} obey an equation of state. 
This equation of state connect them together. Generally, it is difficult to attribute precise vari- 
ables to a thermodynamic system under observation, and sometimes the choice seems subjective. 
From De Donder's work on irreversible chemical reactions, it is known that a supplementary ther- 
modynamic variable is necessary to completely describe, thermodynamically, a system which has 



been brought out of equilibrium 
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2l| . This variable, denoted ^, was called the degree of ad- 
vance of the chemical reaction. Subsequently, it has been generalized by Prigogine, Defay and 
van Rysselberghe to any irreversible process. In this generalization, C, characterizes the advance of 
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the process. This app roach has been used until the emergence of the thermodynamic concept of 



internal variables 
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23|. While in chemical reactions, C, can be well-defined, in other irreversible 



processes, its nature is more difficult to precise. In the glass community, C, has been called the 



degree of order or the order parameter 



46|. Following De Bonder, the variation of entropy of a 



system undergoing a thermodynamic transformation is given by: 

dS = deS + diS (2) 

where d^S is the reversible exchange of entropy between the system and its surrounding (exchange 
of heat, work or matter) and diS is the positive generation of entropy produced within the system 
itself when it is relaxing towards its equilibrium state. The entropy production term and the order 
parameter are connected by the following equality: 

d.S = ^>Q (3) 

T is the temperature of the system and A the affinity of the process. A is the intensive conjugated 
thermodynamic variable associated to the order parameter ^. In the {p, T} state ensemble, the 
affinity is defined as: 

-(f). 

where G is the Gibbs free energy of the system. The affinity is the thermodynamic force of the 
process exactly like gradients of intensive variables for other processes. From the definition of the 
state function G {G = H — TS), the Berthelot-De Bonder formula is obtained: 

-Kf).-(f). 

where H is the enthalpy. A system out of thermodynamic equilibrium is consequently character- 
ized by the variables {p, T, ^} or {p,T, A}and not only by {p,T} variables. The order parameter 
obviously has an equilibrium value for any values of the couple {p,T}. In this case, this equilib- 
rium value ^(;g{p,T), defines the condition of thermodynamic equilibrium for the system with the 
necessary and sufficient condition that A = at any instant. But in the case of the glass transition, 
if the system is (for a certain time) brought out of thermodynamic equilibrium, the three latter 
variables must completely define the instantaneous state of the system. The system is not anymore 
constrained to remain in the {p, T} equilibrium plane. Instead, it leaves this plane to follow the 
value of the affinity and order parameter, which become now time-dependent. ^ and A depend 
of time and of the instantaneous time course of p or T. For example, if the system is suddenly 
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brought out of thermodynamic equihbrium by means of quick pressure or temperature jumps then 
after the arrest of these perturbations, the variables ^ and A continue to evolve until they reach 
their equilibrium values {C,eq for the order parameter and for the affinity). This time evolution 
is expressed in the physical ageing process of glasses. A glass is a system for which the order 
parameter has not any possibilities to reach such equilibrium values, because experimentally the 
system is brought in a range where the relaxation time becomes too high and ^ is frozen. During 
the vitrification process, the heat capacity and other thermodynamic coefficients become the sum 
of two separated contributions due to time-scale separation. The first contribution is concerned 
by the rapid degrees of freedom (phonon's bath) as compared to the speed of variation of the 
temperature. The internal exchanges of energy which follow a temperature perturbation can occur 
between these rapid modes. In other words, the modes of vibration are thermalized. They define 
the temperature T of the system. The other heat capacity contribution involves the slow internal 
degrees of freedom as compared to the speed of variation of the temperature (it is assumed that 
these are connected to the configurational changes inside the system). They are driven by the 
variable ^, and this heat capacity contribution is called the configurational contribution j47|: 

Upon this decomposition, the glass transition is the thermodynamic transformation represented 
by the gradual freezing of these configurational modes when the pressure or the temperature vary. 
That is to say that the second contribution of the right hand side of Eq. ([6]) is progressively 
lost during the time course of p or T. The adjective "progressively" or "gradual" means that this 
loss of configurational contribution occurs in given pressure or temperature intervals, and not 



only at a definite temperature or pressure (like in a phase transition) 19|, |2^. If the principal 
characteristic of the glass transition is a continuous freezing of the configurational degrees of 
freedom, then the variable ^ which is associated to these modes is progressively slow down until 
total arrest. In this case, the velocity of ^ is completely neglected over the temperature time 
rate {{j^)^ = and thus C^"^^ = 0). Consequently, a glassy state is a state defined by a heat 
capacity (idem for the other coefficients) at constant order parameter = {^) ^ where only 
the phonon contribution is measured. Contrary, the heat capacity of a liquid or a supercooled 
liquid contains, in addition to the latter, the order parameter contribution which in this case is 
maximum ($C_{p}"{conf}=Max$). In this case, the internal exchange of energy that follows 
the temperature variation allows the system to adapt its structure to the new temperature. This 
maximal contribution for each thermodynamic coefficient can be exactly calculated as a function 
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of the thermodynamic derivatives such as $(see following sections). We will see that these maximal 
contributions, that are also named $ are those that in fact are computed in the expression of the 
classical PD ratio. 

B. Working assumptions: single order parameter, single relaxation time and linearity 

Here we briefly asses the important assumptions used in the following: 

• Firstly, we suppose that one single order parameter is sufficient to explain most of the basic 
features of the glass transition. An important consequence of such an assumption is that the 
configurational contributions of the three thermodynamic coefficients are all driven by the 
same order parameter C,'- 



Let us also indicate that in the following, others relaxational effects are not taken into 
account. Secondary /3-relaxation processes or anharmonicity effects are not tackled because 
they are not accessible by Cp, Kt and ap measurements. 

• The second assumption requests that the order parameter relaxes with one well-defined 
relaxation time r, although it is known that a distribution of relaxation times is in fact implied 
for a more adequate description of the glass transition. We can, however, suppose that r 
is the average upon all the Tj implied in the distribution. With the first assumption, this 
implies that enthalpy and volume evolve along time following the same r after temperature 
or pressure disturbances. 

• Eventually, we suppose that the glass transition can be well depicted using linear physics 
of irreversible processes. This is the case if temperature or pressure time rates imposed on 
the system are so slow that the system never departs too far from equilibrium. This latter 
assumption leads up to important mathematical consequences: 

1. -affinity and the time rate of the order parameter are proportional: 




(7) 




(8) 




(9) 



8 



where L is the so-called Onsager coefficient. It is sometimes discussed in the literature that 
such a relation can hold even for high departures from equilibrium. For example, the Fourier 
law remains valid even upon high temperature gradients. Prigogine and co-workers have 
experimentally validate such relations in the case of particular chemical reactions where the 



assumption comes back to neglect A as regard to RT (R is the ideal gas constant) j25[. Now, 
discussing the validity of such linear relation for the case of the glass transition is a problem 
beyond the scope of this article. This has been well investigated by Moller and al. and 



Gutzow and al. 



24 



211 



2. -the isobaric and isothermal heat and volume advancements of the transformation (derivat- 
ives of the ffist order with respect to C, of the enthalpy and the volume) have always constant 
values, i.e their values at equilibrium: 




dV 
dH 



eq 



P,T 

eq \ ' 



3. -the second order derivative of the Gibbs free energy with respect to ^ has also its value at 
equilibrium: 



4. -the order parameter relaxes towards its equilibrium value following a simple first order linear 
differential equation: 

- Ceg) ..ON 
~ [16) 

at T 

The equilibrium value of the order parameter C,eq is exclusively and instantaneously driven 
by the temperature or the pressure. 

5. -The last consequence is certainly the most difficult to justify. It states that close to equi- 
librium the affinity is negligible as compared to the isothermal heat of the transformation: 

In Appendix A, we show briefly that this last assumption comes to neglect the heat of friction 
due to the entropy production term as compared to the heat released or supplied due to the 
advance of the order parameter at a given temperature. In other words, this means that 
the amount of heat involved in the structural change is greater than the heat produced by 



9 



Al 
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internal dissipation. This implies also that the Active temperature Tj can be, to a first 
order, close to the temperature T of the system in the linear regime (but bear in mind that 
dTf 7^ dT to a first order). 

these assumptions have been discussed in classical non-equilibrium thermodynamics books 



-|30|. 



III. CONFIGURATIONAL PRIGOGINE-DEFAY RATIO 

A. Derivation of the configurational Prigogine-Defay ratio for general out of equilibrium 

transformations 

A thermodynamic transformation occurs out of equilibrium if the system of interest is, over 
a certain time, departed from its state of equilibrium. This means that the order parameter ^ 
is deviated from its equilibrium value ^eq whereas it evolves along time following the equation 
(fT3|) . The flux of the order parameter is driven by the thermodynamic force of the process which 
corresponds to its conjugated affinity. Since the affinity is a state function, it can be differentiated 
with respect to the three state variables, namely {p, T, ^}. This yields the total differential of 
affinity |l|: 

fdA\ ,^ fdA\ , fdA\ ,^ 
dA= [ — ] dT+ [ — ] dp+ [ — ] d^ (15) 



\dTj,,^ \dPj^^^ V5e/,,T 
The first and second partial derivatives involved in this equation can be transformed using Maxwell 
relations at constant pressure and temperature respectively. With the help of Eq. ([5]) this yields: 



Using the definition of affinity, the third partial derivative is simply (with an opposite sign) equal 
to the second order derivative with respect to ^ of the Gibbs free energy at constant pressure and 
temperature: 

the total differential of the affinity (Eq. (IT^ ). it is now possible to extract the isobaric temperature 



The thermodynamic stability of the system requires that must be strictly positive. From 
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derivative of the order parameter and the isothermal pressure derivative of the order parameter 
involved in the three configurational contributions (Eq. (jS]) and (^): 



dT 



dp y rp 



rp(dA\ 

KdTJp 



(19) 



dV 
95 



T,p 



dp J rp 



W)p,T 



(20) 



These two temperature and pressure derivatives reveal the affinity production due to temperature 
perturbation under isobaric conditions and due to pressure perturbation under isothermal con- 
ditions respectively, which are the signatures of irreversible processes taking place during such 
perturbations. Now, from these three configurational contributions (Eq. ([7]), ([8]) and ([9])): 

(21) 

(22) 
(23) 







Akt = 


conf 


Aap = 


^ conf 
OL 



a configurational Prigogine-Defay ratio can be derived: 

1 ACpA^T _ 1 (^)p,T (f ): 



■j-j- conf 



VT fAa) 



T (dv\ (ii) 



(24) 



Here the sign "A" signifies a difference between thermodynamic coefficients corresponding to the 
effective state of the system measured during the glass transition and those of the corresponding 
glassy state (iso-order-parameter state). With Eq. f[T^ and fl20]) the CPD ratio can be expressed 
as the following: 



■Qcon/ 









(f) 






p,T ^'^^'P 



(25) 



This is the general expression of the CPD ratio during a non-equilibrium glass transition. Some 
comments are necessary at this point. Firstly, the CPD ratio is defined over the whole transform- 
ation interval (temperature and pressure intervals) within which the glass transition takes place, 
and not only at a predetermined temperature such as Tg or a characteristic pressure within the 
glass pressure range. This remark has already been outlined by Schmelzer and Gutzow 19|, |24| . 
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Secondly, due to the time dependence of affinity (and its derivatives), n™"-^ can be time dependent. 
However, time dependent quantities in Eq. fl25]) can vary in such a way that the ratio remains 
fixed and time independent. The notion of time dependent PD ratio has already been envisaged by 
Havlicek 
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32|. Let us, for closing this part, indicate that the assumptions previously mentioned 



have not been used for the moment. Now, from the general expression of the CPD ratio, different 
cases can be envisaged. 



1. Ideal glass transition: usual Prigogine-Defay ratio for an equilibrium thermodynamic transformation 

A transformation which takes place at thermodynamic equilibrium is a transformation for which 
the control parameters (pressure or temperature) vary so slowly that the system has enough time 
to restore its equilibrium state. In other words, when C,eq{p,T) moves, ^ follows it instantaneously, 
or in practises, ^ moves more rapidly than the variations of the control parameters. Consequently, 
it is impossible to observe an equilibrium glass transition. If it would be the case, the liquid 
must remain a super-cooled liquid during all the pressure or temperature courses until the liquid 
becomes a solid. The heat capacity of the super-cooled liquid can be extrapolated to a characteristic 
temperature for which the liquid and glassy states are the same. This temperature is the Kauzman 
temperature, Tk, for which the relaxation time is infinite. To the best of our knowledge, nobody has 
never been able to perform experiments in order to approach sufficiently close to this characteristic 
temperature. As of today, it is unknown whether or not at this temperature a phase transition 
occurs, or whether there a discontinuity in the derivatives of the Gibbs free energy. During a real 
glass transition, eventually, the relaxation time of the order parameter becomes so high that it is 
impossible for the system to adapt its structure to the new temperature and pressure, respectively, 
in the experimental time scale. Thus, a glass transition occurs outside thermodynamic equilibrium. 
Anyway, let us suppose for instance that such an ideal glass transition is possible at the limit of 
infinitely slow temperature and pressure rates: dT/dt — and dp/dt — > 0. In that case, the three 
configurational contributions ACp, Akt and Aop are at their maximum and always equal to the 
difference between the equilibrium liquid and the glass at each temperature and pressure (the glass 
being the state of matter obtained for infinitely fast temperature and pressure rates: dT/dt oo 
and dp/dt — )■ oo). On Figure 1, a schematic diagram shows how the heat capacity behaves during 
different thermodynamic transformations, i.e during an equilibrium transformation (liquid), during 
a perfect non-equilibrium transition (ideal glass) and during a real glass transition (viscous glass). 
At this point, one can ask whether measuring such differences in the thermodynamic coefficients 
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does make sense for the study of a real glass transition. Indeed, even in the presence of a real glass 
transition (for which dT/dt takes a defined value different from zero), if the difference between 
the super-cooled liquid and the glass is taken at Tg by extrapolation of the heat capacity of the 
liquid and the glass at this temperature, this comes back to take ACp'^ for determining the PD 
ratio. Anyway, we can calculate these maximal configurational contributions by means of the total 
differential of affinity in putting A = and dA = 0. In this situation, temperature and pressure 
derivatives of the equilibrium value of the order parameter become (Eq. [T9|20|) : 

dO'' 



dH 
95 



dT 



dp J rp 



eq 



(26) 



9S J 



eq 



9^2 



eq 



(27) 



The equilibrium configurational contributions (Eq. [3, [HI El) are as follows: 

2 



9h\ 



eq 



P,T 



T 



(28) 



P,T 







[(t; 


't,p_ 



V 



dV_ 



( 9!G^ 



eq 



P,T 



eq 



di) 



eq 



P,T 



VT 



92 G 



eq 
P,T 



(29) 



(30) 



We can compute these three equilibrium contributions in the formula of the PD ratio (or simply in 
putting A = and dA = in Eq. ( l25l) ) in order to obtain a PD ratio, in its usual meaning, which 
is strictly equal to unity. The conclusion is that, under the single order parameter assumption, 
during an ideal glass transition the PD ratio is equal to unity like in the case of a second order 
phase transition. During this equilibrium transformation, the PD ratio is measured over a given 
temperature or pressure interval and not only at a critical temperature or pressure. Now a pertinent 
question could be to know if, under the single order parameter assumption, the non-unity of the 
PD ratio could be due to the real out of equilibrium nature of the glass transition. 
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2. Isoaffine glass transition 

The isoaffine transformation is a very particular case of non-equilibrium transformation. The 
fact that this type of transformation can well-describe a glass transition is, in our point of view, 
questionable. Indeed, it must be very improbable that all the partial derivatives in the equation 
(fT5|) are connected in such a manner that dA equals zero at any time during the transformation 
while A is different from zero. Moreover, since we know that the glass transition is a transformation 
starting from an equilibrium state (super-cooled liquid) for which A = to a glassy state for which 
A ^ 0, then the isoaffine glass transition is, in fact, impossible. Nevertheless, it seems that this 
has been the way which was initially followed by Prigogine and Defay, and that this way was also 
recently reconsidered by Schmelzer and Gutzow [l^. A critical lecture of the famous paper of 
Davies and Jones (Ref. jsl) reveals that most of the mathematical thermodynamic derivations are 
also made under the assumption of iso-affinity. Indeed, in making dA = in equation (125|) then 
the expression of the PD ratio derived by Schmelzer and Gutzow is found: 

(f ) 



95 



Under the assumption of linearity (close to equilibrium transformation) then for the isoaffine glass 
transition, n*^""-^ is very close to unity because in this case the inequality A « (^^^ holds (see 
Eq. [m see also more precision on this assumption in Appendix A). 



3. A three dimensional diagram of the glass transition 

The state of a material out of thermodynamic equilibrium is given by three independent ther- 
modynamic variables [p, T, A). Thus the state of a glassy material can be represented in a diagram 
with three axes representing the values of p, T and A as in Figure 2. The order parameter is a 
function of these three variables alone: ^ = ^{p, T, A) while at equilibrium, the equilibrium order 
parameter = C,eq{p,T,A = 0) is fixed by the condition A = and is a function of p,T only. 
Upon cooling at constant pressure p = Po but with different temperature rates 71 and 72, the 
system is driven out of equilibrium by the temperature variation rate and enters into the slowly 
relaxing states at different temperatures Tq and Tf, > Tq. Upon further cooling, the system enters 
a true glassy state with a fixed order parameter and These evolutions are represented as 
trajectories (a) and (b) on the Figure 2. Although glass 1 and glass 2, represented by the point (a) 
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and (b), are at the same temperature and pressure, it is clear that they are not in the same states 
since they have different affinities. Nevertheless if these two glassy materials were still within the 
glass transition interval, they would relax toward the A = and eventually reach identical equi- 
librium states. The fictive temperature concept is also a measure of the structure of the system 
like ^ (cf. following paragraphs). According to the definitions given by Davies and Jones jSj, the 
fictive temperature Tf of the materials can be deduced in the linear regime by the intercept of 
the tangent of the trajectories in the {pq,T,A) plane with the {pq,T,A = 0) axis. In the glass 
transition interval, the fictive temperature is varying and follows the control temperature although 
with some delay, its value at a particular point within the glass transition range of trajectory (b) 
has been represented on the diagram of figure 2. As soon as the system is in a glassy state but still 
in the linear regime, the fictive temperature is constant and noted Tj and necessary the trajectory 
in (po,T,A) plane becomes rectilinear. The final fictive temperatures (Tj) have been marked on 
the diagram for trajectories (a) and (b). Next we perform a cooling of the material under pres- 
sure pi at a rate 72 as represented by trajectory (c). We explicitly supposed in the diagram that 
dTg/dp > (where Tg is the glass temperature). In the case of a true first order transition, this 
would be equivalent to supposing that the molar volume of the glass is smaller than the molar 
volume of the liquid for a positive melting latent heat. In other words, this means that under pres- 
sure for a given temperature rate, the system is driven out of equilibrium at higher temperature 



than under atmospheric pressure 



33I-35I. This can be re-phrased in terms of fictive temperature: 



for a given cooling rate, the fictive temperature of a sample cooled under high pressure will be 
higher than the fictive temperature of a sample cooled under atmospheric pressure. This has been 



recently observed experimentally 



33N35|. The final fictive temperature has also been marked for 



this glass. Finally we perform an experiment in which the system is maintained in equilibrium 
{A = 0) at constant temperature Tq. The pressure is then increased at a given rate yielding the 
trajectory (d) pictured on figure 2. If the pressure variation drives the system out of thermody- 
namic equilibrium, then above a vitrifying pressure pg, the system enters first a glass transition 
range in pressure and if the pressure is further increased, the system becomes totally glassy and 
increasing the pressure does not change significantly the value of the order parameter ^. So, if the 
cooling rate of trajectory (c) and the increasing pressure rate of trajectory (d) are well chosen, it 
is possible that the two trajectories intercept as pictured by point (d) in the diagram. A fictive 
pressure of the glassy material (d) can be defined in the same manner as the fictive temperature. 
From this diagram, it can be seen that the system can be set into an identical glassy state by using 
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either temperature variation at constant pressure or pressure variation at constant temperature. 
This naturaUy provides a hnk and equivalence between the Active pressure and Active temperat- 
ure concepts and the order parameter approach introduced by De Bonder. Figure 3 provides an 
experimental illustration of above arguments. It shows the dependence of a structural parameter, 
the relative amount of tetrahedrally coordinate boron, BO4, in a borosilicate glass as a function 
of fictive pressure and fictive temperature 33|. Here, [-BO4] may be taken as related to the order 



parameter, and different {T,p) paths may lead to different structural states of the final glass. In 
the following, a formalism is developed that will help to translate future physical experiments in 
the framework of the CPD ratio. 



B. Expression of the Configurational Prigogine-Defay ratio as a function of the isobaric 
heat capacity and isothermal compressibility coefficient 



From now, the assumptions described in paragraph 2.2 will be used in order to simplify the 
general expression 025 p of the CPD ratio. This will allow to connect its value to measurable 
quantities in a different, novel fashion. Firstly, let us re- formulate expression (l25l) in the form 
of a ratio of two adimensional quantities that involve only thermal and mechanical variables, 
respectively: 

jjconf 



n 



calo 



where we have for the two terms of the ratio: 



n 



calo 



•4+(f),,-^(#), 



dv\ 



dV 



dp J J, 



Using the assumptions of paragraph 2.2, we obtain: 



n, 



calo 



^« )p,T 



dvY" 



(32) 



(33) 



(34) 



(35) 



16 



It could be easily demonstrated that the two previous quantities can be expressed as a function 
of the heat capacity only for Hcaio snad as a function of the compressibility only for Hmeca (see 
Appendix B): 

"^"^'^ ~ ACV ^ ^ ^^^^ 

TT /\ K.rj-i /\ K.rp 

T^meca ~ = (38) 

ACp"^ is the equilibrium contribution as classically involved in the usual PD ratio although ACp is 
the difference between the effectively measured heat capacity during the glass transition (viscous 
state) and the glassy state. The same arguments hold for the measured isothermal compressib- 
ility. Thus, the CPD ratio now depends only on the isobaric heat capacity and the isothermal 
compressibility (the isobaric thermal expansion coefficient does not now play a role anymore): 

neon/ ^^p" ^I^T ( . 

Let us remark that Ilcaiois what is known in the glass literature as the inverse of the normalised 
heat capacity which is equal to unity only when the system is in the liquid state. In letting time to 
the material to adjust its configuration at each pressure and temperature, respectively, we observe 
that H-caio and Hmeca tend directly toward unity. The question is whether this attainment holds 
on such a manner that n^°"-^ remains equal to one at each instant. This could be envisaged in 
developing the latter expression of II'^"'^'^ using the fact that ACp and Akt involved temperature 
rate and pressure rate, respectively: 

' p,T dT/dt ^\dO p,T dT/dt ^^^^ 

A., = -lf^l d_^^_^(dVY d_^ ^^^^ 



y J Tp dp/dt y J Tp dp/dt 
With these previous equations, in assuming that near equilibrium, a temperature variation under 
isobaric condition induces the same order parameter disequilibrium as does a pressure variation 
under isothermal condition, with the help of (Eq. (l28l) and fl29|) ) we obtain for n^°"'^ the following 
value: 

av\ 
^^)t,p 



^conf . 1 \9i)p,T dT/dt 

TfavY" dp/dt ^ ' 



17 



This provides a condition for the unity of the CPD ratio under relaxation involving mechanical 
and thermal variables and temperature and pressure time rates: 



eg 



dV 

dT/dt _ V^Ot,p .... 
dp/dt fdHY' ^ ^ 

This is a prerequisite for experiments in order to obtain a CPD ratio equal to unity. The ex- 
perimental challenge is to be able to follow this requirement for being at the same point in the 
three dimensional diagrams of figures 2 and 3 after temperature and (or) pressure perturbations. 
Noteworthy, as mentioned before, if enthalpy and volume relax towards equilibrium with different 
relaxation times for equivalent thermodynamic perturbations, the ratio of the two relaxation times 
{th and Ty) must appear in the previous equation. In other words, if thermal and mechanical 
variables are driven by different order parameters, then experimentally, we must adapt the rate of 
variations of T in one experiment and the rate of variation of P in another one in order to fulfil 
the unity condition for the CPD ratio. Wondraczek and collea gues have performed experiments in 



order to probe such P,T influences on the structure of glasses 



33 



C. Expression of the configurational Prigogine-Defay ratio as a function of the Active 

temperature and the Active pressure 



Tool has originally defined a concept, called the fictive temperature, adapted to a phenomenolo- 



gical description of calorimetric experiments on glasses |36l. l37l|. In short, this quantity represents 
the temperature of the slow relaxing modes during a glass transition. Moynihan and co-workers 
have used the fictive temperature concept for a description of the glass transition and ageing pro- 



cesses, and have connected it to heat capacity data |38H40[. By means of a systematic conformity 
between the "two temperatures thermodynamics" developed by Nieuwenhuizen and the order para- 
meter approach of De Donder, we have recently demonstrated that the fictive temperature of Tool 
can in fact represent a real thermodynamic temperature for systems undergoing irreversible pro- 



cesses 



41] . For instance, the fictive temperature of Tool is equivalent to the effective temperature 



of Nieuwenhuizen. Davies and Jones have however already established such a direct link between 
the fictive temperature and the order parameter {s]. In the same way, they also have defined a 
fictive pressure. Close to equilibrium, the affinity can be developed as a function of the fictive 
temperature departure from the classical one (under isobaric condition) and also as the function 
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of the fictive pressure departure from the pressure (under isothermal condition): 

A-{T-Tf)[^—^^^ (44) 

^~(p-,,)(-j^^^ (45) 

With equations f|T6|) and (fT7|) and neglecting the affinity as compared to the heat of the transform- 
ation, we obtain: 

Using linear relation (fTOj) this yields to a new expression for JY^""^^ from Eq. (l42l) : 

neon/ ^ (P - P/ ) ^^/^^ ^48) 

(T - T^) dp/dt ^ ^ 

Using the well-know differential equation of Tool and a similar one for the fictive pressure in the 
linear regime we have: 

ATr (T - TA 

(49) 



dTf _ (T - Tf) 
dt T 
dpf _ {p~Pf) 



dt T 

When T or p are stopped, these two linear differential equations give an exponential decay for the 
fictive temperature and the fictive pressure which are linked to enthalpy and volume relaxation. As 
we are placed in the assumption of one single order parameter (r is the same in the two equations), 
this leads to a simple expression for the CPD ratio: 

It is straightforward to remark that a result obtained a long time ago by Moynihan, connecting 
the normalised heat capacity with the temperature derivative of the fictive temperature has been 
38, 4o||: 



re-found 



U^aio ~ ^ = ^ (52) 

Here ACp"^ is taken at T while in the original derivation of Moynihan ACp"^ is taken at Tf in the 
normalised heat capacity expression. The two ways are equivalent upon the linear assumption 
because T ~ T/- (see Appendix A for details). The conclusion of this paragraph is that the CPD 
ratio can involve only fictive temperature and fictive pressure along with the classical temperature 
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and pressure. Notwithstanding its simplicity, this sort of connection can be a good way to estabhsh 
the hnk between this theoretical approach and physical experiments. Relation (15T]) can however 
indicate to us how, firstly the CPD ratio must be equal to unity or not, and secondly, how the 
classical PD ratio, which is the long-term average of n™"'^, can also be equal to unity or not. This 
is briefly discussed now. 

IV. ON THE UNITY OF THE CONFIGURATIONAL PRIGOGINE-DEFAY RATIO 

At infinitely long time-scale, i.e for infinitely slow cooling rate or pressure variations, the system 
has enough time to restore its internal equilibrium. It behaves like a liquid. The fictive temperature 
and pressure vary exactly like their corresponding temperature and pressure driving variables. Eq. 
S5T\i states that n^°"-^ is equal to unity. But as already mentioned, under these circumstances, no 
glass transition does occur. In real situations, even for very slow variations of the driving variables, 
there exists an instant at which the relaxation time of the system r(T, p) will reach a value for 
which the system cannot anymore restore its long range configuration. From this instant, the 
experimenter begins to observe a fall in the measured thermodynamic coefficients {Cp, etc.). This 
fall is not sudden. It occurs over a definite interval (temperature or pressure). During this collapse, 
the system looses its configurational entropy which is observable by a smooth jump in the heat 
capacity. This being, there is also relaxational effect which becomes small due to the importance 
of the relaxation time. There is however an irreversible generation of entropy. This entropy gain 
due to friction during slow internal configurational readjustments is generally small. In the linear 
regime, the entropy production term is indeed negligible (see Appendix A). This has already been 
discussed by Davies and Jones jsj). From this situation, the fundamental question arises as to 
whether the departure of the fictive temperature from the equilibrium temperature of the system 
reaches a value in such a fashion that the fictive pressure attains a value different from the pressure 
of the system such that the Eq. (15T]) leads to unity for 11™"-^. Indeed, from Eq. (HHj) . considering 
that only one single order parameter is implied, ]J'^°'^f is unity if the scanning temperature or 
pressure rates (under isobaric and isothermal conditions respectively) are in proportion with the 
departure of the fictive temperature and pressure respectively from their driving variables: 

dT/dt ^ (r - r^) ^^g^ 

dp/dt (p — Pf) 

Here, the requirement obtained in paragraph 3.2 (Eq. fj43|) ) is found as a function of the fictive 
temperature and the fictive pressure directly, which are measurable quantities. The discussion 
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made after the Eq. (H3|) is still valuable. The fact that the CPD ratio should be equal to unity only if 
;his type of statement is fulfilled is an argument in favour of researchers (see for example McKenna 
2!) who argue that the departure from unity for IT remarked in the literature is probably due to a 
erroneous consideration of the different histories in the measurement of the three thermodynamic 
coefficients Cp, and a^. However, the statement that "11 different from unity" is due to the 
presence of more than one single order parameter is also true as we have shown that in this 
paper. But, to our point of view, if we suppose that "11 different from unity" is probably due to 
a bad consideration of the history of a glass sample in the measurement of the thermodynamic 
coefficients, then "11 different from unity" is not a definitive proof in favour of a multi-order 
parameters description of the glass transition. 



V. CONCLUSION 



In this paper, we have shown that, when the configurational parts of the three thermodynamic 
coefficients involved in the expression of the PD ratio are considered, a general expression can be 
derived for this ratio. We have called this general expression, the configurational PD ratio. The 
classical PD ratio usually discussed in the literature, and experimentally measured by extrapol- 
ation of the properties of the super-cooled liquid and the glass at T^, is a limiting case of this 
expression for infinitely high time scales of observation (or identically for infinitely slow temper- 
ature and pressure time variations), i.e when only thermodynamic equilibrium is considered. We 
have questioned the pertinence of the classical expression of the PD ratio for the characterization 
of the glass transition since it is known that the glass transition is by essence an out of equilibrium 
transformation. Since when classical PD ratios are measured, the obtained values are generally 
different from unity, we have discussed that it could be due to the out of equilibrium nature of the 
glass transition. Indeed, in a real glass transition the configurational parts of the three coefficients 
Cp, and Op vanish gradually upon temperature decrease or by varying the pressure. Some 
classical and some more recent studies have tackled this issue in discussing the PD ratio in terms 
of difference of the thermodynamic coefficients between a state at constant affinity and a state at 
constant order parameter. We have also criticised this point of view, because during the out of 
equilibrium phase transformation there is a permanent creation of affinity and the isoaffine case is 
not typical for the glass transition. If the out of equilibrium character of the glass transition must 
be correctly taken into account then we must construct a PD ratio with thermodynamic coefficients 
effectively recorded during the glass transition, i.e by taking into account a difference between an 
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out of equilibrium state during the transition (measured Cp, kt and ap) and a state at constant 
order parameter (the glass: Cp^^, Kt,£, and ttp,^). When taking into account such time dependent 
variations of the three thermodynamic coefficients, under the single order parameter assumption 
and under the assumption of identical thermodynamic histories, the unity for the CPD ratio is 
always confirmed. Our conclusion at this level joins that which can be found in the literature on 
this topics, i.e, that the non unity of the classical PD ratio is either due to the presence of several 
order parameters for a complete description of the glass transition, or because in most of the cases 
ACp, Akj- and Aap are measured on sample for which thermodynamic history have not seriously 
be taken into account For the moment the debate is not closed. Under linearity conditions, 
we have also given an expression for the CPD ratio where the number of relevant coefficients have 
been reduced from three to two, and where only normalised isobaric heat capacity and normalised 
isothermal compressibility play a role. From this last expression the CPD ratio can be expressed 
as a function of not only the temperature and the pressure but also as a function of the Active 
temperature and the Active pressure. All these properties are experimentally accessible, and a 
future challenge will be to correlate the present formalism to experimental situations. 
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Equilibrium 

transformation: Cp |jqLijd(/5-^0) 




Figure 1: On this schematic graph Cp as a function of temperature upon coohng (/3 = dT/dt) is shown 
in different situations. An equihbrium transformation occurs when the coohng ramp is so slow that 
equihbrium is maintained ah the time, i.e the measured heat capacity is that of a super-cooled liquid. A 
completely non-equilibrium transformation occurs when the cooling ramp is infinitely fast. The system be- 
comes instantaneously frozen, i.e the measured heat capacity is that of a glass. A real transformation takes 
place when a liquid is cooled following a determined scanning rate where it undergoes a measurable jump 
in the heat capacity. The difference ACp'^ which is computed in the classical PD ratio is taken between 
the hypothetically equilibrium and completely non-equilibrium state taken at the glass temperature Tg of 
the real glass transition. 
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A 



A 




Figure 2: Diagram in the {p, T, A) space representing the state of an out of equihbrium system. Different 
hypothetical trajectories yielding glasses are represented: (a) p = po; ^ = 7i; (b) P = Po'i ^ = 72 > 7i; 
(c) p = pi > poj ^ = 7ij (d) ^ = ^0) ^ 7^ 0- The corresponding glass states are marked by the letters 
(a), (b), (c) and (d). The equilibrium state from which the glasses are formed is marked by the letter 
(e).are connected by equilibrium (reversible) trajectories in the (p,T,A = 0) plane. 
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Figure 3: {Pj, Tj, [S04]}tri-dimensional graph showing the dependence of a structural parameter, the 
relative amount of tetrahedrally coordinate boron, BO4, in a borosilicate glass as a function of fictive 
pressure and fictive temperature 



33| . The concentration of boron [i?04] is the structural order parameter 



in these experiments. 
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Appendix A: 



One assumes that sufficiently close to equilibrium the affinity obeys the inequality: 

^«f^y' (54) 



Dividing by T and multiplying by d^/dt the two members yields: 

M « i^ii^ (55) 

T dt T ^ ^ 

The left-hand-side is the so-called entropy production term and the right-hand-side the configur- 
ational entropy rate: 

diS dS^'^'f 

dt dt ^ ^ 

The configurational entropy rate is defined thanks to the configurational heat capacity such as 

follows: 

^Sconf _ X f 

~ T ^^^^ 

Thus neglecting the affinity as compared to the heat of reaction at constant temperature and 
pressure is equivalent to neglect the entropy production as compared to the configurational entropy 
change. 

Davies and Jones have shown that close to equilibrium the affinity can be written simply as 
follows jsl: 

^~'^-^''y., '''' 

This comes simply from a Taylor series expansion with an arrest to a first order in the departure of 
the fictive temperature to the temperature. Using Maxwell relation and the Berthelot-De Bonder 
formula (Eq. [S]), this last equation becomes: 



(59) 



Let us first observe that if the inequality (iMj) is fulfilled then the affinity can be written: 



Secondly, if we develop further the previous equation fl59l) . without using the inequality flM|) . we 
arrive easily to the following equation: 
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In comparing the two previous equations, we conclude that neglecting the affinity before the heat 
of reaction in the linear regime is also equivalent to state that the Active temperature remains 
close to the real temperature of the system to a first order: 

Tf-T (62) 

Let us be aware, such as shown in the paper, that it is not true for the first order derivative of 
these two temperatures even to a first order: 

dTf ^ dT (63) 

Knowing this fact, we can start from the expression of the PD ratio recently derived by Schmelzer 



and Gutzow for isoaffine case 



19|: 



And using the equality (|6T!) we obtain easily that: 

n^Le - ^ (65) 

which, as demonstrated, in the linear regime is close to the unity. Such an expression for the PD 



ratio has already been obtained by Ottinger 



by an other way 



42|. If we start mostly from the 



equality (l60l) . that we have derived recently 41|, we obtain for the isoaffine PD ratio the other 
expression: 

T 

n™=Le ~ 2T-Tf ^^^^ 
which is also very close to the unity in the linear regime. 

A last remark can be made: in using the Eq. (l6T]) in place of Eq. fl60|) in the Eq. fH2|) then, a 
more exact expression is found for the expression (15T]) which becomes now: 

^conf _ TidpidT_ ^ dpf dlriT ^^^^ 
T dp dTf dp dlriTf 

Upon the assumption used, this is equivalent to the expression ( I5T]) . 



29 



Appendix B: 

Firstly, let us rewrite the exact expressions of Hcaio and n„ 



n 



dH 



calo 



A + 



dH 



^ \dT)i 



9y\ 



dV_ 



T,p 



dp J J, 



(6^ 



(69) 



For the equation (168|) . in multiplying the numerator and denominator by the ratio 
(^^) ' obtain the expression: 



n 



calo 







2 

/T 





















(70) 



where we recognize for the numerator an expression close to the maximum contribution of the 
configurational heat capacity AC^^ (see Eq. 1281) . Upon the assumptions used, where it is admitted 



that 



eq 



^] ' and 



P,T 



d^G 



eq 



P,T 



then it is equivalent to AC^'^. For the 



denominator, we recognize ACp expressed in function of the affinity and the temperature derivative 
of the affinity (see Ref. 43|). Indeed, in using the total differential of the affinity under isobaric 
condition, we can easily replaced [d^/dT)^ in the classical expression of the heat capacity: 

di 



C -('^ 



dH\ fdH 



dT 



(71) 



in order to obtain for ACp 
equilibrium there is: 



Cp^^, which is the desired expression. In conclusion, close to 



n 



calo 



AOT 



(72) 



Starting from the equation (l69l) . the same reasoning can be done. Let us start from the definition 
of the compressibility coefficient: 

~ [dPj^ " V Vd^)T,p [dpJT 
In using the total differential of the affinity under isothermal condition we have obtained: 



(73) 
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which included in the preceding equahty yields to: 

2 



dV 



+ 



dV 
95 



dp J 



T 



V ( ^ 



(75) 



P,T 



We observe that such expression appears directly if we multiply the numerator and denominator 
of the equation ([69]) by {%) ^ /V (^0^ ^ , and in this case the maximal contribution of the 
configurational compressibility coefficient (in the linear regime) appears directly: 

2 



A 4^ 





r" 




T,p 



V ( ^ 



eg 
P,T 



(76) 



Consequently: 



^meca ~ 



(77) 



